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The nonlinear magneto-plasmonics [1—3] aims to utilize plasmonic excitations to control the mech¬ 
anisms and taylor the efficiencies of the non-linear light frequency conversion at the nanoscale. We 
investigate the mechanisms of magnetic second harmonic generation in hybrid gold-cobalt-silver mul¬ 
tilayer structures, which support propagating surface plasmon polaritons at both fundamental and 
second harmonic frequencies. Using magneto-optical spectroscopy in Kretschmann geometry, we show 
that the huge magneto-optical modulation of the second harmonic intensity is dominated by the ex¬ 
citation of surface plasmon polaritons at the second harmonic frequency, as shown by tuning the 
optical wavelength over the spectral region of strong plasmonic dispersion. Our proof-of-principle ex¬ 
periment highlights bright prospects of nonlinear magneto-plasmonics and contributes to the general 
understanding of the nonlinear optics of magnetic surfaces and interfaces. 

Rapid development of plasmonics has facilitated an outstanding progress in understanding, designing and control¬ 
ling the optical response of metallic nanostructures, including that in nonlinear-optical domain [3-0. The optical 
second harmonic generation (SHG) from solid interfaces represents a well-known experimental technique with numer¬ 
ous applications in physics, chemistry and biology. Being a natural tool to enhance the light-matter interaction, the 
plasmon-assisted localization of the electric field in small volumes proved to be very effective at elevating the effi¬ 
ciency of the nonlinear-optical processes. Beyond the physics of plasmonic nanoantennas requiring the fabrication of 
sophisticated metallic nanoobjects [8. 9], the Kretschmann geometry for the excitation of surface plasmon polaritons 
(SPPs) in a thin metallic layer on a dielectric prism plays a special role due to its robustness and simplicity [loj. A 
huge plasmonic enhancement of the electric field caused by the phase-matched excitation of SPPs under the resonant 
coupling conditions has been shown to boost the SHG intensity EMI- The most recent experiments have evidenced 
that not only the SPP-resonance at a fundamental frequency uo but also the excitation of a second harmonic SPP at 
the doubled frequency 2c o plays an important role in the process of the nonlinear-optical conversion EMI. However, 
due to their dispersive nature a simultaneous excitation and coupling between the SPPs at fundamental and second 
harmonic frequencies has not been reported so far. 

The main idea of this study was to design a nanoscale optical experiment where we could excite, detect and ma¬ 
nipulate the nonlinear interactions between the SPPs at the fundamental and SHG frequencies. To achieve this goa l 
we have adopted the concepts of linear magneto-plasmonics in hybrid metal-ferromagnet multilayer structures p, [l5| 
together with the most recent ideas in nonlinear magneto-plasmonics EHEIE3. Our experimental geometry is 
sketched in Fig. la. A thin gold/cobalt/silver trilayer structure was grown on a glass substrate by means of the mag¬ 
netron sputtering. A 5 nm-thin magneto-optically active layer of ferromagnetic cobalt was protected from oxidation 
by a 3 nm-thin layer of gold. A 25 nm-thick silver layer acted as main constituent in this hybrid plasmonic nanos¬ 
tructure, which was excited by collimated femtosecond laser pulses through the glass prism (Kretschmann geometry). 
The reflected SHG intensity, as well as the linear reflectivity signal at both fundamental uj and double frequency 2c o 
was recorded as a function of the incidence angle 0 for the two opposite in-plane directions of magnetization (+M 
and — M in Fig. la) of cobalt in the transverse Kerr geometry. 

Understanding the properties of the SPP dispersion [l8| in our structure (see Fig. lb and Supplementary Information 
for details) is one of the key points for nonlinear magneto-plasmonics. Due to inevitable dispersion, the SPP excitations 
at fundamental and double frequencies in Kretschmann geometry occur at slightly different angles, suggesting that 
simultaneous phase-matched excitation of both SPPs is impossible. However, nonlinear-optical considerations account 
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FIG. 1: Surface plasmons in Au/Co/Ag multilayer structures, a Sketch of the experimental geometry showing the 
spatial distribution of the square of the normal projection of the electric field at the angles of incidence 6 corresponding to the 
SPP excitations at the fundamental uj and double 2uj frequencies. Insets in the coloured circles show the experimental angular 
dependencies of the linear-optical reflectivity and SHG intensities for the two opposite directions of magnetisation of the cobalt 
layer, b Dispersion of the SPP in the Au/Co/Ag trilayer under study. Black dash and green solid lines represent the photons 
in vacuum and in glass, respectively. Thick dot purple line is the calculated SPP dispersion, its linewidth is shown with the 
purple background area. The inset illustrates the possibility of a simultaneous excitation of the SPPs at both frequencies u 
and 2u. 


for the following nonlinear phase-matching condition 

2 k 0 (u)n(cj) sin 9 n \ = fc spp (2w), (1) 

between the excitation source at the silver-glass interface characterized by the in-plane component of the k- vector 
A# sin 0 and the second harmonic SPP at the gold-air interface with the A;-vector k% p p . Being just one of several possible 
SPP frequency conversion pathways Bin, this phase-matching condition plays a crucial role in our interpretation 
as it determines the resonant SPP-induced enhancement of the nonlinear susceptibility • 

The experimental results are summarized in Fig. 2. The linear optical reflectivity in Kretschmann configuration 
using the fundamental (920 nm) and second harmonic wavelengths (460 nm, frequency doubled in a BBO-crystal) 
are displayed in Fig. 2a. A narrow plasmonic reflectivity dip at 920 nm observed at 42.5 degrees contrasts with a 
much broader minimum at 460 nm at 6 ~ 47 degrees, in agreement with much higher SPP-losses at 2a;. In terms of 
linear magneto-plasmonics these two resonances behave similarly: the relative magnetisation-induced variations of the 
reflectivity, A R/R neatly follow the angular differential reflectivity dR/dO , suggesting that altering the magnetization 
direction induces a small shift of the reflectivity spectrum due to the magnetic contribution to the SPP wavevector 
110 . 

Whereas the linear reflectivity at fundamental frequency uj shows only small modulations of the order of 1%, in 
line with Ref. [6], the SHG angular spectra at 2uj display drastic changes, as shown in Fig. 2b for four different 
values of the fundamental wavelength 760, 800, 820 and 920 nm. Red and blue background areas show the SPP 
dispersion (see Fig. lb) at the frequencies uj and 2a;, respectively, recalculated as a function of incidence angle 0. 
Over the SHG wavelengths range (380-460 nm) the SPP shows a strongly dispersive behavior resulting in a shift 
of the resonant angle and decrease of losses (see the linewidth). Fundamental SPP remains nearly dispersionless 
and demonstrates no noticeable change of its linewidth. Contrary to previously reported results on a gold film 0 , 
the angular positions of the SPP resonances for the fundamental and SHG frequencies in our multilayer structure 
correspond to the pronounced minima in the SHG intensity. A strong dependence of SHG intensity on magnetization 
direction is quantified by the magnetic SHG (mSHG) contrast p: 

= h U (+M) - I 2 u (-M) 

^ + l2u( — M) 

shown in Fig. 2d. Surprisingly, the largest magnetisation-induced modulation of the SHG intensity is accompanied 
by the SPP excitation at the SHG frequency and not the fundamental one. Angular dependence of mSHG contrast 
re-plotted in Fig. 2e without an offset displays the largest reported values of modulation reaching 33 % around 
9 = 44 degrees. For the shortest wavelength (760 nm) the mSHG maximum is only about 20%, since in this case the 
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FIG. 2: Angular spectra of the linear and nonlinear magnetooptical response in presence of the SPPs. a,b 

Linear-optical reflectivity R and magnetization-induced reflectivity variations A R/R for the excitation with the 920 and 460 
nm wavelength. c,d SHG intensity and (d) mSHG contrast p angular spectra plotted with an offset. Red and blue background 
areas illustrate the SPP dispersion in the experimental spectral range (760-920 nm for fundamental, 380-460 for the SHG, 
respectively.) e Angular spectra of the mSHG magnetic contrast p for the four fundamental wavelengths, no offset is introduced. 
Inset: Illustration of the magnetisation-induced changes of the total SHG intensity, huj{±M) oc \E 2u} (±M)\ 2 = \E 2u} d= E^\ 2 . 


SPP damping at the SHG frequency becomes so large that the system approaches the region with the non-propagating 
(propagation length A SPP < A 2 J) SPPs. This observation links up our measurements to the most recent results by 
Zheng et al. [2], who reported similar values of the mSHG contrast on a 10 nm-thin iron film on glass, a structure not 
supporting SPPs at the SHG frequency. This fact, along with the dispersive shift of mSHG maximum reinforces our 
conclusion that the 33 % large mSHG contrast (which is equivalent to the increase of the SHG intensity by a factor 
of 2 upon magnetization reversal) is dominated by the SPP resonance at the SHG frequency. 

In order to understand these experimental observations we have used the approach proposed by Palomba and 
Novotny 0 , who explained the complex angular dependence of the SHG intensity from a thin gold film on glass by 
an interference of two contributions coming from the gold-air and gold-glass interfaces. In our case (Fig. 3a) the silver- 
sapphire interface acts as a source of the nonmagnetic SHG , and the upper part consisting of Au and Co layers 
is assumed to generate the electric field E^ containing both magnetic and non-magnetic x^ contributions. 

Thus the total SHG intensity hu is described by: 

I 2u oc I El“ + E 2u ± E 2 Z I = |xf } : E X E X + X $, 2) : E 2 E 2 ± X < 2m) : E 2 E 2 1 2 , (3) 

where complex tensor components Xi and xF ±xf m ^ represent the effective optical nonlinearities at both interfaces. 
The angular dependences of the electric fields E\ and E 2 at the fundamental frequency a;, which are driving the SHG 
process at both interfaces, are calculated by finite difference time domain (FDTD) method for the 860 nm excitation 
wavelength (see Fig. 3c-d). 

Based on these calculations, the modulation of SHG intensity in the vicinity of fundamental SPP resonance can 
be explained by destructive interference of two contributions Ef^ and E^ 0- For our multilayer structure the 
amplitude of the SHG light generated at the upper interface and transmitted through the prism at any angle 6 turns 
out to be smaller as compared to the other one: |F? 2u; | < \Ef UJ \. This is why the strong increase of E 2 in Fig. 3c-d at 
the top interface explains the SHG intensity minimum around the fundamental SPP resonance at 0 = 42.5 degrees. 

The same explanation holds for the SHG minimum around 6 = 45 degrees albeit for a different reason. We note 
that the nonlinear susceptibility x® is known to acquire a SPP-induced resonant contribution, x® = Xnr + Xrel 
00 with: 


XresW ^ 


1 

0 — 0 n \ + iT 


( 4 ) 


Owing to this resonant contribution, the SHG field enhanced at the top interface destructively interferes with the one 
generated at the bottom interfaces, which explains the experimentally observed SHG intensity minimum. In order 
to understand the angular spectrum of this resonant contribution we approximated it with the Lorentzian resonance 
line. The width T and the angular position 6 n \ of this resonance line was chosen according to the calculations shown 
in Fig. lb for the nonlinear SPP excitation at frequency 2a;, or 430 nm wavelength. This resonant shape successfully 
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FIG. 3: Schematics of the SPP-assisted mSHG generation in a multilayer structure, a Illustration of the origin 
of the inequality of the SHG intensities for the opposite direction of the external magnetic field H leading to the non-zero 
magnetic contrast p. b SHG sources at the two interfaces, bottom (1, glass/Ag) and top (2, air/Au/Co/Si). c-d False colour 
spatial distributions of the square of the electric field projections \E Z \ 2 and \E X \ 2 at the 860 nm fundamental wavelength. 
Dashed white lines show the sample borders. Solid white lines represent the corresponding \E\ 2 angular distributions at the 
top and bottom interfaces, e-f SHG intensity and mSHG contrast angular spectra (open circles and squares) for the 860 nm 
fundamental wavelength together with the fit lines based on the Eq. ©• Dashed lines are the angular profiles corresponding 
to the resonant SPP excitation at the frequency uj and nonlinear SPP excitation uj —>• 2uj (see Fig. 5S in the Supplementary 
Information). 


accounts for the large mSHG modulation at the angles corresponding to the nonlinear SPP excitation quantified by 

Eq. (1). Based on Eq. (3) we were able to fit the experimental angular spectra using the resonant part of the nonlinear 
(2) 

susceptibility Xres only. 

Leaving out the details of our fitting of Eq. (3) for the Supplementary Information, where the choice of predominant 
Xcomponents (from the six non-zero ones, three magnetic and three non-magnetic, at each interface) also justified, 
we would like to claim to have obtained a good quantitative agreement between the theory and experiment both for 
angular dependence of the SHG intensity (Fig. 3e) and mSHG contrast (Fig. 3f). Note that the mismatch between 
the experimental data and the fit curves occurs largely around the fundamental SPP resonance, where the resonant 
Xcontribution used in our simulations is hardly playing any role. 

Based on our numerical calculations, we can conclude that the large value (33%) of mSHG contrast is dominated 
by the properties of SPP resonance at SHG frequency. Thus, upon reversing the external magnetic field the SHG 
output can be doubled, which opens up a new strategy for the design of nonlinear magneto-photonic devices operating 
at the nanoscale. The frequency dependence over the dispersive SPP spectral range clearly demonstrated that the 
maximum value of 33% mSHG contrast remains large whereas its angular width decreases as the SHG resonance gets 
narrower and shifts towards the fundamental resonance. 

In the given spectral range one would expect to reach further enhancement of mSHG contrast by systematically 
varying the individual thicknesses in this trilayer structure. The future route for investigations will imply moving 
towards the experiments in the telecom frequency range, which apart from the obvious technological importance would 
allow to extend the magneto-optical investigations including the magnetization-induced third harmonic generation. 
Moreover, expanding into the far-IR or THz spectral range, where the SPP dispersion is lower, would allow to explore 
the phase-matched magneto-plasmonic coupling between the SPPs at the fundamental, second, third and higher-order 
harmonic frequencies. With the largest reported plasmon-assisted value of mSHG contrast of 33% we get closer to 
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the dream of the nonlinear magneto-photonic devices based on novel physical principles. 
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Supplementary information 


Surface plasmon dispersion in Au/Co/Ag trilayers 


In order to obtain the surface plasmon polariton (SPP) dispersion relation 


0 


^spp(^) — &spp ~f~ ^spp — ^o(^spp ^spp) — 


1 + £eff 


( 5 ) 


which is characterized by frequency-dependent SPP wavevector k spp (uj) or refractive index n spp = n' pp + m' pp we 
applied the effective medium theory. It allowed us to approximate the effective dielectric constant of the Au/Co/Ag 
interface 0: 


1 

£ e ff = t- / e(z)e 

Cskin JO 


z /s s ^dz 5 


(6) 


where the dielectric function e(z) depends on the normal coordinate 2 inside the multilayer structure. The spectral 
dependence of the dielectric function for gold, cobalt and silver is taken from Ref. [3]. The total thickness of two 
upper metal layers (3 nm Au and 5 nm Co) is significantly smaller than the penetration depth of the electric field 
S s kin ~12 nm. 

In order to explain the results of the graphical representation in Fig. lb we would like to quantitatively compare 
the SPP dispersion relations at Ag/air and Ag/5nm Co/3nm Au/air interfaces. Figure IS shows the wavelength 
dependence of the real (n' pp ) and imaginary (n' pp ) parts of the SPP refractive index as well as the SPP skin depth 
and propagation length L spp . The presence of a highly absorbing cobalt layer in our multilayer structure leads to a 
significant reduction of the SPP propagation length L spp = l/(2kon” pp ) over the entire spectral range, as compared 
to the Ag-air interface. In contrast to this behavior, the skin depth £ s ki n and the real part of SPP refractive index 
n' pp exhibit substantial differences only in the vicinity of the SPP resonance corresponding to an optical wavelength 
in vacuum of 350 nm. Moreover, in the entire spectral range of our SHG measurements (380-460 nm for SHG and 
760-920 nm for fundamental wavelength), the skin depths for both structures appear to be nearly identical. This 
allows us to calculate £ e ff in Eq. © using the skin depth for the Ag-air interface. 
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Figure IS: Spectra dependence of complex refractive index n spp = n' pp + in " pp , skin depth £ s kin and propagation 
distance L spp for the SPPs excited at the air-Ag and air-Au/Co/Ag interfaces. 

Now we turn to the discussion of Fig. 1 in the main manuscript. Due to substantial losses the (linear) resonance 
condition in Kretschmann configuration 


k x {0,u) = fc' pp (w) (7) 

with k x (6,uj) = ko(uj)n(uj) sin 6 appears to be fulfilled within a finite range of incident angles 0 (or tangential compo¬ 
nents k x (Q,uj) of optical wavevectors in the prism). 

To account for this angular broadening, we have applied a simple phenomenological model, where the complex 
Lorentzian-shaped function f(k x ,uj) determines the efficiency of the SPP excitation: 

f f i \ _ ^sppM 

/( “ ,w) “ (fc'ppH-M + ^'ppH’ ( } 


or, equivalently 


/(#><*>) 


_ fc "ppH _ 

( fc s PP M - ko{u)n{u) sin6>) + ifc" pp (w) ‘ 


( 9 ) 


The shaded areas in Fig. lb and Fig. 2c,d in the main manuscript are limited by two contours of this angular 
distribution at the half-maximum of this Lorentzian \f(k x , uj)\ 2 or |/(0, cj)| 2 . 
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Linear reflectivity in presence of SPPs 


Angular dependence of the linear reflectivity in Kretschmann geometry was measured for eight different wavelengths, 
namely 760, 800, 860, 920 nm at fundamental laser frequencies and 380, 400, 430 and 460 nm with its second harmonic 
light (generated in a BBO-crystal before the sample), see Fig. 2S(a,b). Maximum values of the reflectivity are 
normalized to unity because these supplementary measurements were not calibrated to get absolute values. Angular 
spectra of the magneto-optical Kerr effect measured at the four fundamental frequencies (Fig. 2S,c) demonstrate great 
similarity to each other due to low SPP-dispersion in this spectral region. 

Despite the fact that the curves measured in the blue (2a;) range also look somewhat similar, the change of the 
resonant dip when decreasing the wavelength is still noticeable. Figure 2S(c,d) illustrates angular spectra of the 
transverse Kerr effect at the same wavelengths. Note the characteristic asymmetric shape of the spectra, present for 
both red and blue sets of wavelengths. It is worth mentioning that the transverse Kerr effect at blue wavelengths is 
several times smaller than at the red ones, which obviously correlates with the quality factor of the SPP resonance. 
Furthermore, the magnitude of the magnetoplasmonic modulation of the reflectivity decreases when the wavelength 
approaches 380 nm. This observation can be explained by comparing the propagation length of a SPP L spp with 
the light wavelength A (see Fig. IS). It is seen that L spp decreases together with the fundamental wavelength, thus 
reducing the quality of the SPP excitation and also the Co magnetization-induced contribution to the SPP fc-vector. 
As the latter is directly responsible for the SPP-assisted magnetoplasmonic reflectivity modulation, the transverse 
magneto-optical Kerr effect also decreases, albeit a very broad dip in the reflectivity is still seen. 
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Figure 2S: (a-b) Angular spectra of the SPP-modulated reflectivity R for the four red (a) and blue (b) wavelengths, 
(c-d) Angular spectra of the magneto-optical transverse Kerr reflectivity modulation dR/R for the same wavelengths. 


Nonlinear-optical SPP excitation 

A conventional way of exciting a SPP in the linear optics requires the linear phase-matching condition © to be 
fulfilled [l]« This linear phase-matching condition describes an SPP excitation at a frequency uj. However, this phase¬ 
matching condition appears to be modified, if an SPP is generated via the nonlinear frequency mixing, i.e. surface 
second harmonic generation (SHG) at the frequency 2 uj [H, 5]. A straightforward analogy with SHG generation in a 
nonlinear crystal directly leads to the following phase-matching condition 

k x ( uj) + k x (uj) = 2ko(uj)n(uj) sin# = k' spp (2u) 


( 10 ) 
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for two incident photons generate generate an SPP at doubled frequency. As in case of linear resonance this phase- 
matching condition is broadened by SPP-dissipation 

f nl (k x , 2 W ) = (fc , pp(2w) _ 2M W) ) + ifc" pp (2 W ) • (U) 

Due to very high plasmonic losses the acceptance angle for phase-matched SPP generation at 2uj is much larger 
as compared to the linear phase-matching at the fundamental frequency. Note that this nonlinear phase-matching 
condition slightly differs from the linear phase-matching condition for a collimated beam at frequency 2uj\ 


f(k x , 2w) 


_ fcspp(M _ 

( fc spp( 2w ) - k x (2uj)) + ifc"pp(2w) ' 


( 12 ) 


Comparing the expressions for 2 k x {uS) = 2kon(uj) sin 6 and k x {2uj) = 2kon(2uj) sin# we recognize that the angles 
for the phase-matched SHG generation and linear excitation in Kretschmann configuration are different only due to 
dispersion in a glass prism: n{uj) ^ n(2uj). Knowing this small angular shift of about 0.7 degree between f n \(6,2uj) 
and f(6 , 2uj) for 860 nm fundamental wavelength, we can correctly calculate the SPP-assisted field enhancement at 2c u 
from the linear reflectivity in Kretschmann configuration. As we will see below, this analysis appears to be particularly 
useful to estimate the angular dependence of the resonant part of the x^ -tensor. 


Fitting experimental data 

General considerations 

Trying to identify the leading tensor contributions in second harmonic generation from a isotropic surface is a 
delicate issue of playing with three possible contributions (consider p-polarised fundamental and SHG beams): Xzzz , 

Xzxx and Xxzx @. In the chosen Cartesian coordinate frame z-axis is perpendicular to the surface plane and x-axis 
lies in the incidence plane. 

There are multiple considerations regarding the appropriate choice of the nonlinear susceptibility components, if 
their values are a priori not known. The trajectory of free electrons moving in the direction perpendicular to the 
interface (^-direction) is expected to exhibit the largest nonlinearities due to the structural changes across the interface 
(a large built-in electric field at the interface explaining the physical origin of the work function/contact potential) 
0. Therefore the xizz component is expected to dominate in a free-carrier approximation, when the energy of SHG 
photon is smaller than the energy separation E g between the Fermi level and the d-band. This is indeed the case 
for silver, as E g = 3.8 eV exceeds the largest photon energy used in our experiments (3.3 eV). As described in the 
main Manuscript, there are two sources of SHG in our model, which are located at the Ag/glass and Ag/Co/Au/air 
interfaces. Given that all our attempts to get reasonable fits with the components other than x\ izz Ag/glass 
interface failed, we fixed this zzz-component in further simulations described below. We attribute the ability to 
identify the dominant tensor component to the aforementioned fact that both relevant frequencies co and 2 uj are below 
the d-band edge in silver. 

The above conclusions largely rely on the ongoing debates regarding the choice of the most important components 
of the nonlinear susceptibility x ^ [ffl. Naraoka et al. [T(| demonstrated that Xzzz yields a better agreement 
with the experimentally observed SHG output enhancement. The same conclusion is implicitly confirmed by other 
authors 7, [Tl| . However, several authors [13, !!l3j consider instead the Xxzx component, following the original paper 
by Simon et al. 0, which nonetheless claims the similar behaviour of all three components in the vicinity of the 
SPP resonance. Note the rigorous, although somewhat less instructive approach used by Pavlov et al. 0 , where all 
components were taken into account, and thus the number of the fit parameters increased dramatically. We attribute 
these discrepancies to a simple observation that it is generally difficult to isolate the dominant contribution to x^ 
tensor if either the fundamental or SHG-frequencies overlap with interband transitions. 
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Figure 3S: (a) Angular dependence of electric fields at the gold-air interface. Whereas the tangential component E x 
is continuous across the interface, the normal component E z experiences a large jump, (b) Dependence of electric 
field on the normal coordinate z at # = 42.5° corresponding to the SPP excitation at the frequency uj. 


However, different contributions P 2uj = XijkEjEk to the non-linear polarization also depend on the components of 
the local electric field Ej^. The boundary conditions in Maxwell equations dictate the continuity of the tangential 
component of the electric field E x and the normal component of displacement field D z = s(z)E z across the interfaces. 
In our plasmonic case characterized by \eau\ ^ 1 the boundary conditions suggest E z >> E x in the dielectrics (air and 
glass) and E z <C E x in Ag, Au and Co. This statement is confirmed by Finite Difference Time Domain calculations of 
the electric fields for our structure leading to \E X | 2 ~ 25\E Z | 2 for electric fields inside the metal and thus explaining why 
we have obtained better fits considering the X^ x ^x contributions. We would like to note for the sake of completeness 
that on the air side the relation is opposite: \E X \ 2 — 0.04 \E Z \ 2 . This suggests that the nonlinear-optical response 
of our plasmonic structures covered with a nonlinear medium (such as, for example, a monolayer of molecules or 
semiconductor quantum dots) will be probably dominated by different contributions to the optical nonlinearity. 


Fitting SHG intensity 

As a first step in the fitting procedure we decided to try different x ^ components at the top interface, in attempt 
to figure out which one is the most important for the observed modulations of the SHG-signal. In order to make 
the fitting procedure more transparent we have limited ourselves to the following simplified versions of Eq. (3) in the 
main Manuscript: 


J zzz « \X ( llzz E lz s[ne + X < 2 l i ) zzz E lz s[ne \ 2 
!xzx « \X^zzz E l,z sin $ + X^]czx E 2,x E 2,z COS #| 2 

! 2 zxx « IXl 2 L- E l 2 ,z sin6 ' + X2 2 Lx^h sin6, | 2 > 


(13) 

(14) 

(15) 


where Fd, E 2 are the calculated electric fields at the bottom and top interfaces, respectively, magnetization-induced 
^( 2 ,m) components are neglected, and X^Izzj Xj2xzx an d X^lxx are th e complex non-magnetic components of the 
tensor. Additional factors sin# and cos# originate from the geometrical projection of the normal (z) and tangential 
(x) components of the nonlinear polarisation on the propagation direction of the SHG-beam in the glass prism (see 
Fig. 4S). 
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Figure 4S: Geometrical interpretation of SHG in Au/Co/Ag trilayer structure on a glass prism. We assume that SHG 
occurs at the bottom (1) and top (2) interfaces. SHG electric fields in the glass prism depend on SHG polarization 
(normal, linearly polarized along the z-direction or tangential, oscillating along the x-direction), which results in a 
additional geometrical factors sin# and cos#. 

These equations can be further simplified in order to minimize the number of the fit parameters: 


t2uj 

1 zzz 

OC 

\ E 1,Z sin 0 + r zzz El z sin#| 2 

(16) 

t2uj 

1 xzx 

OC 

l-Ei,* sin 9 + r xzx E 2 , x E2,z cos 9 2 

(17) 

t2uj 

1 zxx 

OC 

\ E 1,Z sin 0 + r zxx E l tX sin 9\ 2 

(18) 


Without the loss of generality we have used complex coefficients r zzz , r xzx and r zxx which represent the ratios of the 

( 0 \ 

respective yh^-components to Xi, ZZ z ^ the experimental data. The results of the best fitting shown in Fig. 4S (a) 
were obtained using r zzz = 0.04 • exp(il.057r), r xzx = 0.19 • exp(il.67r) and r zxx = 0.9 • exp(i0.057r), respectively. Three 
dashed contours illustrate the Lorentzians given by expressions (18111112ft . The nonlinear phase-matching condition 
shifts the position of SPP-resonance for SHG-generation to larger angles as compared to the linear SPP resonance at 
2cj. 
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Figure 5S: Angular dependence of SHG intensity for 860 nm optical pump together with best fitting results based on 
nonresonant Xnr model (a, see Eq. ([TS1 14) and resonant Xres($) model (b, see Eq. ([2T1 19H. Three dashed Lorenzians 
at the top of the figure aim to illustrate the angular position and width of the SPP-resonances at fundamental and 
SHG frequencies. 


Two main features can be recognized in Figure 5S(a). First, in the vicinity of the fundamental resonance (0 ~ 42.5°) 
all three expressions adequately reproduce the experimental data and giv e almost identical fits. Therefore, our fitting 
attempts support previously mentioned difficulties and ambiguities [14J in the identification of the dominant x ^ 
components at arbitrary metal surfaces. 

Second, the behavior of our fit curves around 2uj SPP resonance (0 ~ 45 — 46°) is notably different from the 
experiment. This failure to reproduce the experimentally observed SHG intensity minimum stems from the fact that 
in the above models we have assumed that x ^-tensor components do not depend on the angle 6. This assumption 
appears to be fundamentally wrong when x ^ contains resonant (SPP-induced) contributions at the SHG frequency 
2 cj, thus highlighting one of the most important features in our modeling of the SPP-assisted SHG. 

The usual way to take into account SPP-resonance at SHG-frequency 2 cj is based on the representation of x ^-tensor 
components as a sum of the non-resonant and resonant contributions [6|, 7]: 

X (2) = +Xres( 6 ')' ( 19 ) 

In our case the resonant contribution originates from SPP-excitation at SHG-frequency 2 uj. Whereas the attempts to 
fit the experimental data shown in Fig. 5S(a) were based on the non-resonant contribution x ^ — Xn? (which does 
not depend on #), here we focus on the resonant part Xres(#)- We assume that the angular dependence of the resonant 
contribution is as well governed by a Lorentzian shape , which can be obtained from Eq. (j9j) by linearizing sin# in 
the vicinity of 0 n \: 



r r ( res ) v ( 2 ) 

e-e nl + ir 


(20) 
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Here T = k' s ' pp (2u:)/(ko(uj)n(uj) cos^ni) = 2.8°, and 0 n \ = 45.8° is given by Eq. (6). The expressions for SHG intensity 


with resonant contributions 


y (2) 

Ares 


thus read: 


t2u 

1 zzz 

oc 

Kz 

t2uj 

1 xzx 

OC 

Kz 

T 2(jJ 

1 zxx 

(X 

Kz 


sin # + f nl {Q)El z sin #| 2 
sin # + f nl (0)E 2 , x E2,z cos #| 2 

sin# + f n i(9)E% tX sin6»| 2 . 


( 21 ) 

( 22 ) 

(23) 


The results of the best fitting in Fig. 5S(b) indicate that all three equations with = 0.06 • exp(*0.37r), rxzx = 

0.25 • exp(z0.87r) and rzxx = 1-35 • exp(zl.257r), respectively, provide a better approximation to the experimental 
data. However, the I^ xx demonstrates the closest proximity to experimental data thus corroborating our qualitative 
conclusion that the contributions containing larger local electric fields \E X \ \E Z \ should dominate. 


Fitting mSHG and magnetic contrast 




Figure 6S: Angular dependence of SHG intensity for 860 nm optical pump together with best fitting results based on 
resonant Xrll(0) model (see Eq. ([241 -22)) for two opposite orientations of the magnetization M. 

With the results of the two previous sections at hand, we are now ready to fit our experimental SHG intensities and 
mSHG contrast to the Eq. (3) in the main Manuscript. We note that in the final expression for the SHG intensity 
we already know the field distributions E XiZ , which account for the SPP resonance at the fundamental frequency c j, 

as well as the resonant contribution of the nonlinear susceptibility Xres($) associated with the SPP resonance at the 
double frequency 2uj. 

In order to stay in line with systematic fitting procedures we have considered all three magnetic nonlinear sus¬ 
ceptibility components x£zT\ X^zzx^ and x£xx l) E3 leading to the following expressions for the SHG intensities 
/(it M): 


Ilt m) ex | E 2 hz sin 9 + r^Ue)[El x sin# ± cos#]| 2 (24) 

Izzx m) \ E lz s[nd + r zxx fm(S)[El x smd ±r^E 2 , z E 2 , x sme ]\ 2 (25) 

I 2 ^ ex \E 2 hz sin# + r( r x e x s )/ n i(#)[ J F 2 2 x sin# ± cos#]| 2 . (26) 

Figure 6S shows the final results of our fitting in a direct comparison suggesting that Ixxx provides the best fitting 
results, in agreement with our interpretation of dominating terms proportional to \E x (uj)\ 2 . The angular spectra of 
the magnetic SHG contrast are presented in Fig. 7S. 
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Angle 0 (degrees) 

Figure 7S: Angular dependence of mSHG contrast for 860 nm optical pump together with best fitting results based 
on resonant model with three different magnetic components of %res-tensor(see Eq. (JT51 -9)). All fits accurately 
reproduce the the behavior of the mSHG contrast in the vicinity of the SHG-resonance and, notably, the sign change 
at 0 n \ = 45.8 degrees. 

The fitting results in Fig. 5S and 6S were obtained with a fixed, previously deduced value r^xx = 1.35exp(il.257r) 
and using the normalized magnetic components /X^cx as fit parameters: r^l = 0.02 • exp(—i0.47r), 

r^x = 0.07 • exp(i0.057r) and = 0.25 • exp(i0.57r). We note that among these values r^xx is in a good agreement 
with the results of theoretical calculations for Fe predicting a value |r m | ~ 0.25 [18]. 

The behavior of mSHG contrast in Fig. 7S is accurately reproduced in the vicinity of the SHG-resonance at 
0 n \ = 45.8° but still shows significant deviations at the fundamental SPP-resonance around 42.5°. We attribute 
these differences to the fact that our simplified model disregards the non-resonant contributions to the x^ -tensor, 
the inclusion of which, however, would inevitably stir up troubles related to the multiparameter fitting. Systematic 
experimental studies of similar multilayer structures with different layer thicknesses should provide more insight about 
the details of the magnetic SHG and mSHG contrast behaviour at the fundamental SPP-resonance. 
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